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An effective field theory technique that combines the standard nuclear physics approach and chiral 
perturbation theory is applied to the hen process, "^He + n — > *He+7. For the initial and final nuclear 
states, high-precision wave functions are generated via the variational Monte Carlo method using 
the Argonne v\4, potential and Urbana VIII trinucleon interactions, while the relevant transition 
operators are calculated up to 0{Q'^) in HBxPT. The imposition of the renormalization condition 
that the magnetic moments of "^He and be reproduced allows us to carry out a parameter-free 
calculation of the hen cross section. The result, a — (60±3±1) fib, is in reasonable agreement with 
the experimental values, (54 ± 6) fib and (55 ± 3) fib. This agreement demonstrates the validity of 
the calculational method previously used for estimating the reaction rate of the solar hep process. 



One of the important advantages of using effective 
field theory (EFT) in nuclear physics is that it leads to 
well-controlled predictions for nuclear observables over 
which the conventional model approaches have little or 
no control. One such case is the hep process '^He -\- p ^ 
"^He-l-z^e+e^, which can be potentially important for solar 
neutrino physics and non-standard physics in the neutron 
sector. In a recent paper, Park et al. ijy have developed 
an EFT approach, to be referred to as "MEEFT" 0, and 
succeeded in making a totally parameter-free calculation 
of the hep cross section. In MEEFT, the relevant tran- 
sition operators are systematically derived using chiral 
perturbation theory, while the initial and final nuclear 
wave functions are generated using the standard nuclear 
physics approach (SNPA). Here SNPA refers to an ap- 
proach in which the wave functions of light nuclei (mass 
number A) are obtained by solving exactly the A-body 
Schrodinger equation with a Hamiltonian that contains 
a high-precision phenomenological potential. With the 
use of MEEFT, Park et al. arrived at a prediction of 
the hep S'-factor with ~ 15% accuracy. Given the fact 
that the predicted values of the hep S'-factor ranged over 
two orders of magnitude in the past, the result obtained 
in P] could be considered a remarkable feat. It is to 
be noted, however, that the ~ 15% uncertainty associ- 
ated with the hep S'-factor (although acceptable on its 
own merit) is strikingly larger than the uncertainty as- 
sociated with the S-factor for the pp fusion process, the 
latter estimated to be less than 0.5% 0. There are two 
plausible causes for such a "large" uncertainty in the hep 
calculation. First, since the calculation in jj] relies en- 
tirely on the information provided by the two-body and 
three-body observables, it may miss part of the dynam- 
ics which is intrinsically four-body in nature. Secondly, 
symmetry mismatch between the initial and final states 
of the hep process, a feature that is absent in the two- 
and three-body systems, induces a strong suppression of 
the otherwise dominant single-particle contribution, re- 
sulting in an enhancement of higher-order corrections, 
perhaps including four-body interactions. 

The purpose of this Letter is to test the reliability of 
MEEFT used to compute the hep process by applying 



the same method to the hen process 



(1) 



at threshold; hen is a four-body process that has close 
similarity to hep in that the pseudo-orthogonality of the 
initial and final wave functions strongly suppresses the 
leading-order one-body contribution. In fact, the degree 
of suppression for hen is such that the exchange current 
"corrections" become dominant terms. Meanwhile, since 
the experimental value of the threshold hen process is 
known with reasonable accuracy: aexp — (54 ± 6) /i6 
and aexp — (55 ±3) fib hen provides a good testing 
ground for the validity of MEEFT. It is to be noted that 
these experimental values have never been explained sat- 
isfactorily before. We remark, in particular, that very 
elaborate SNPA calculations by Carlson et al. and by 
Schiavilla et al. |^ give (7=112 fib and a = 86 fib, respec- 
tively. It has been noted that a is extremely sensitive to 
the ^Re-hn scattering length^a„; the use of the updated 
value, aJfP — 3.278(53) fm 7], makes the calculated of 
a even larger. In view of the high quality of the wave 
functions used in these calculation, it is very probable 
that the problem lies in the many-body currents used. 
This observation makes the use of MEEFT for hen all 
the more interesting. 

Our parameter-free calculation based on MEEFT gives 
for the threshold hen cross section 



cr= (60.1 ±3.2 ±1.0) fib, 



(2) 



which agrees with the above-quoted experimental values. 

The strategy of MEEFT, explained in [J, Q, is close in 
spirit to Weinberg's original scheme 'o'l based on the chi- 
ral expansion of "irreducible terms" . MEEFT has been 
used with great success to describe the highly suppressed 
isoscalar amplitude pertaining to the n -\- p —>■ d -\- ^ pro- 
cess [lol [TTi | and many weak processes (for a recent re- 
view, see Ref. 0). In MEEFT, the relevant current 
operators are derived systematically by applying heavy- 
baryon chiral perturbation theory (HB^PT) up to a spec- 
ified order. The nuclear matrix elements are then eval- 
uated by sandwiching these current operators between 
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the high-precision wave functions obtained in SNPA. In 
EFT, short-distance physics is represented as contact 
counter-terms, and the coefficients of these terms, cahed 
the low-energy constants (LECs) are determined by re- 
quiring that a selected set of experimental data be re- 
produced. This rcnormalization procedure is expected 
to remove to a large extent the model-dependence of 
the short-range behavior of the adopted phenomenolog- 
ical wave functions, and this feature is very important 
for making model-independent predictions. We should 
note here that the strategy employed here is closely re- 
lated to that used for the universal Viow-k in the re- 
cent renormalization-group approach to nuclear interac- 
tions 

The total cross section of the hen reaction reads 
2auj 1 



that the two-body (2B) currents in momentum space are 
valid only up to a certain cutoff A. This implies that 
the currents should be appropriately regulated. This 
is the key point in our approach. This cutoff defines 
the energy/momentum scale of EFT below which re- 
side the chosen explicit degrees of freedom. We use a 
Gaussian regulator in performing the Fourier transfor- 
mation from momentum to coordinate space. All the 
2B currents can be written in the form of Ji2{ri,r2) — 
g-jq (ri+r-2)/2j-^^(^^ _ ^.^y rj,^^ J^^LO corrcctiou is due 

to the so/t-one-pion-exchange. 



Jl2 



where gA = 1.2601, 5*1 



y |(4He|jT(9)|n-H^He;lM}|2, (3) we have also defined f. 
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where a is the fine structure constant, (w, q) is the four- 
momentum of the outgoing photon with |qr| = lo ~ 
20.521 MeV, w„ is the velocity of the incident thermal 
neutron, M is the spin quantum number of the initial 
system, and jV is the electromagnetic (EM) current op- 
erator transverse to the photon direction. As the process 
is isotropic at threshold, we can choose q — z. 

The formalism employed here is essentially the same 
as in the hep process, so we refer to Ref.[l| whenever 
appropriate. We obtain the EM current using HB^PT, 
which contains the nucleons and pions as explicit de- 
grees of freedom, with all the other massive fields in- 
tegrated out. In HB^PT, operators carry a typical 
momentum scale of the process and/or the pion mass, 
which is regarded as small compared to the chiral scale 
A^ 47r/^ ~ TTiAT ~ 1 GeV, where ~ 93 MeV is the 
pion decay constant and m^y — 939 MeV the nucleon 
mass. The leading space-components of the EM currents 
are of order of Q, and hence we refer to terms of order 
Q""*"^ as N"LO. In this paper, we shall limit ourselves 
up to N'^LO or O (Q*) . It is important to take into ac- 
count the kinematical suppression of the time component 
of the nucleon four-momentum and the smallness of the 
photon energy/momentum. For simplicity, we count the 
time component of the nucleon four-momentum and 
q| as order rather than Q. It is important to empha- 
size that, up to N^LO, three-body and four-body current 
operators do not appear. This means that the nuclear 
systems with A — 2,3 and 4 can all be described by the 
same current operators up to this order. 

The leading order one-body (IB) currents are given by 



2miq Amjsi 
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where, p = i(i V — i V) acts only on the wave func- 
tion. Relativistic corrections to the IB currents are 
also included in the calculation. It should be noted 



T,12 — 

T^[ff-CTQ~ io-g] , T0 = Xl T2, (Tq = CTi CTa, 

= ±, X. The explicit forms of the regulated Yukawa 
functions are given in Ref . p| . 

There are no corrections at N^LO. At N'^LO, there oc- 
cur vertex corrections to the one-pion exchange j^^^ , the 
two-pion-exchanges j'^^ , and the contact counter-term 
contributions j^"^. The j^^^ investigated in detail in 
Ref. [13 im reads 
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jtr = X { - |^(c„ + ca) [(4+^ + rf^y-t 

+ (f W + f (->)y2\] - igjCA^'^yJA 
+ iij2NwzTi ■ T2 [(T+y^f^ + (3ff • CT+ - cr+)yj^] }, (6) 

where ~ 0.1021, ca ^ 0.1667, and Nwz ^ 0.02395. 
Although j^"^^ contains the so-called "fixed terms" 
needed to ensure Lorentz covariance, we neglect them 
here since they are expected to be small. The two-pion- 
exchange diagrams, except for the box-diagrams, can be 
found in ref.[l_0|. Adding the box-diagram contributions 
and subtracting the reducible iterated onc-pion-exchange 
part, we obtain 
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The loop functions K' s are defined in Ref. [I[l3. Finally 

there are contact contributions of the form 



•CT 
J12 



2m, 



-q X [34s(cri + 0-2) 



-g.vT^s'YlHr), 



(9) 



3 




FIG. 1: The ^He + n S-wave phase shift in degrees as a 
function of the center-of-mass energy in MeV. The sohd hne 
gives our results obtained with the Woods-Saxon potential, 
while the full dots represent the experimental data obtained 
from an i?- matrix analysis [l^. 



where g4s — mpg4 and g4v — —mp{G^ + jE^'^), nip is 
the proton mass; 34, and E^'^ are the coefficients of 
the counter-terms introduced in A highly note- 

worthy point is that, although there are many possible 
counter-terms, only two linear combinations of them are 
independent for the Ml transition. This reduction of the 
effective number of counter terms is due to the Fermi- 
Dirac statistics; a similar reduction has been noted for 
hep Ul , where only one linear combination of LECs needs 
to be retained. 

As the cutoff A has a physical meaning, its choice is 
not arbitrary. Obviously, A greater than A^ should be 
excluded. Furthermore, A should not be larger than the 
mass of the vector mesons that have been integrated out; 
thus, A < 800 MeV. Meanwhile, since the pion is an 
explicit degree of freedom in our scheme, A should be 
much larger than the pion mass in order to ensure that all 
pertinent low-energy contributions are properly included. 
In the present work, we adopt A = 500, 600 and 800 MeV 
as representative values. 

We use realistic variational Monte-Carlo wave func- 
tions obtained for the Argonne vl4 two-nucleon in- 
teraction and the Urbana VIII trinucleon interac- 
tions; these are the potentials used in the previ- 
ous SNPA calculation 5] . The asymptotic behavior of 
the ■^He + n system has been taken into account by 
means of the Woods-Saxon (WS) potential, Vws[f) = 
Vb (1 + e^''"^"-*/'^'') . We can determine the parame- 
ters in the WS potential so as to reproduce the exper- 
imental data for the scattering length and the low- 
energy phase shifts Q; the results are Vb = 8.63637 
MeV, i?o = 4.14371 fm and do = 0.8 fm. It should be 
pointed out that a small difference a„ can affect the 
hen cross section significantly. For this reason the ex- 
perimental data are directly encoded into theory. The 
resulting phase shifts are shown in figure ^ The precise 
measurement of the ^He -I- n scattering length enables us 
to limit the uncertainty associated with the scattering 
length to < 1% in the matrix elements. The nuclear ma- 
trix elements of the currents are evaluated by taking 1 



TABLE I: The magnetic moment of and ^He (in units 
of Hn), and the hen matrix elements {j)hen (in units of 10~^ 
fm^/^), calculated for A = 600 MeV. The numbers in the 
parentheses represent Monte-Carlo statistical errors. 





m('h) 




{j) hen 


IB 


2.6068(1) 


-1.7839(1) 


-1.76(1) 


2B(l7r + ItvC) 
2B(27r) 
2B(CT) 
2B (total) 


0.2676(3) 
0.0334(1) 
0.0710(3) 
0.3721(1) 


-0.2975(3) 
-0.0334(1) 
-0.0128(3) 
-0.3437(1) 


5.89(3) 
0.90(1) 
0.35(1) 
7.14(3) 


1B+2B 


2.9789 


-2.1276 


5.39(3) 



million samplings in the Metropolis algorithm. 

For each value of A, we adjust the values of 
the two independent LECs (denoted by 54s and 1741,) 
so as to reproduce the experimental values of the 
■^H and ■^He magnetic moments. The results are 
{54s, 544 = {1.079(9), 2.029(6)}, {1.277(12), 0.981(7)} 
and {1.856(22), -0.235(12)} for A = 500, 600 and 800 
MeV, respectively. Once the values of g^sAv ^-re fixed, we 
can predict the Ml transition amplitude for hen without 
any unknown parameters. 

In TableQl we list the contributions of various terms to 
the magnetic moment of and ^He and the hen matrix 
elements {j)hen = |(^He|(j^ - ijy)\n + ^He; 11), calcu- 
lated for A = 600 MeV. The row labelled "2B(l7r-|-l7rC)" 
gives contributions from eas. H5l6|l . "2B(27r)" from eq.jT)), 
and "2B(CT)" from eq.®; "2B(total)" represents the 
sum of all these contributions. As mentioned, the 
IB contribution is strongly suppressed by the pseudo- 
orthogonality of the wave functions, and as a result the 
matrix element is dominated by the (2B). Furthermore, 
since the 2B and IB contributions have opposite signs, 
there is destructive interference between them, which, as 
in the hep case, makes it a more challenging task to ob- 
tain an accurate estimate of the cross section. A highly 
significant feature is that the ratio of the 2B to IB con- 
tribution for hen is about —4, which is much larger in 
magnitude than the corresponding ratio hep, which is 
about —0.6. This large difference can be understood by 
recalling the "chiral filter mechanism" argument 15] ac- 
cording to which the N-'^LO contribution is non- vanishing 
for hen, while the 2B correction in hep only starts at 
N'^LO. It is remarkable that, despite this difference, ex- 
actly the same MEEFT strategy works for both hep and 
hen. 

The cutoff dependence of the matrix elements and the 
corresponding cross section are shown in Table ^ We 
remark that the IB matrix element has no cutoff de- 
pendence. It is noteworthy that the cutoff dependence 
seen in the hen case is quite similar to that seen in the 
hep case. While the individual contributions of the con- 
tact and non-contact terms vary strongly as functions 
of A, their sum shows a greatly reduced A-dependence. 
This can be interpreted as a manifestation (albeit ap- 
proximate) of the renormalization group invariance of 
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TABLE II: The cutoff dependence of the matrix elements 
(in 10~^ fm'''^^) and the corresponding total cross sections for 
A = 500, 600 and 800 MeV. 



A(MeV) 


500 


600 


800 


IB 


-1.76(1) 


-1.76(1) 


-1.76(1) 


2B (non-contact terms) 
2B(contact terms) 
2B total 


5.24(2) 
1.80(1) 
7.04(2) 


6.79(3) 
0.35(1) 
7.14(3) 


8.31(3) 
-0.99(1) 
7.32(6) 


1B+2B 


5.29(2) 


5.39(3) 


5.57(6) 


a{fib) 


56.9(3) 


59.2(5) 


63.2(10) 



physical observables. The smallness of the cutoff depen- 
dence, ^ 3%, seen in Table |n] indirectly indicates that 
our MEEFT scheme allows us to control short-range dy- 
namics to a satisfactory degree. The remaining small 
cutoff dependence may be attributed to the contribu- 
tions of terms ignored in this calculation, e.g.^ the "fixed" 
terms, rt-body currents for n > 2, other higher chiral or- 
der terms, etc. Our final value for the threshold hen cross 
section is cr = (60.1 ± 3.2 ± 1.0) fib, where the first error 
comes from the cutoff dependence and the second from 
the statistical errors. Errors arising from the uncertain- 
ties in the n-^He scattering length are estimated to be 
~ 1% in the matrix element. 

We have already mentioned that the existing SNPA 
calculations for hen cannot explain aexp- A comment is 
in order on this point since, for all the other cases so 
far studied in both SNPA and MEEFT (pp fusion, hep, 
u — d scattering, radiative np-capture etc), the numeri- 
cal results exhibit a close resemblance (except that the 
MEEFT results come with systematic error estimates.) 
As mentioned, one of the most important ingredients of 
MEEFT is a "renormalization" procedure in which the 
relevant unknown LECs are fixed using the experimen- 
tal values of observables in neighboring nuclei. A similar 
procedure has been done in the SNPA calculations of the 
above-quoted cases. However, the existing SNPA calcula- 
tion of hen 1^ 15 , lacks this "renormalization" , and this 



explains why the "existing" SNPA calculation of auen 
disagrees with the experimental data. 

To summarize, we have carried out an MEEFT calcu- 
lation for hen up to N'^LO and provided the first theo- 
retical explanation of the experimental value of the hen 
cross section, cr. Because of the drastic suppression of 
the leading-order one-body term, which makes the sub- 
leading exchange-current a dominant term, and a sig- 
nificant cancellation between these two contributions, an 
unnatural amplification of the cutoff dependence occurs, 
as in the hep case. Despite this accidental "fine-tuning" 
situation, we have shown that MEEFT allows us to make 
a parameter-free estimation of Uhen with accuracy better 
than ~ 10%. The successful appHcation of MEEFT to 
hen renders strong support for the results of our previ- 
ous MEEFT calculation for hep; furthermore, it demon- 
strates the basic soundness of the MEEFT strategy in 
general. The present treatment is open to several im- 
provements such as going to the next order in chiral per- 
turbation, using more accurate SNPA wave functions, a 
more stringent control of mismatch in the chiral counting 
between SNPA and a formally accurate chiral expansion 
that enters in the currents, a better understanding of the 
role that the counter terms play in the renormalization 
group property, etc. It is reasonable, however, to expect 
that the effects of these improvements are essentially ac- 
commodated in the above-quoted error estimate based 
on the cutoff dependence. A robust estimation of the 
hep S-factor has been a long-standing challenging prob- 
lem in nuclear physics . We believe that our MEEFT 
calculation of hep and hen have solved this problem to a 
satisfactory degree. Short of doing a fully formally con- 
sistent EFT calculation, which is at present out of our 
reach, MEEFT seems the best one could do at moment. 
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